Location of files originally: https://github.com/wpetry/USTreeAtlas

library(sp)
library(raster)
library(sf)
## Linking to GEOS 3.13.0, GDAL 3.10.0, PROJ 9.5.1; sf_use_s2() is TRUE
library(maps)
library(tidyverse)
## ── Attaching core tidyverse packages ──────────────────────── tidyverse 2.0.0 ──
## ✔ dplyr     1.1.4     ✔ readr     2.1.5
## ✔ forcats   1.0.0     ✔ stringr   1.5.1
## ✔ ggplot2   3.5.1     ✔ tibble    3.2.1
## ✔ lubridate 1.9.2     ✔ tidyr     1.3.0
## ✔ purrr     1.0.2
## ── Conflicts ────────────────────────────────────────── tidyverse_conflicts() ──
## ✖ tidyr::extract() masks raster::extract()
## ✖ dplyr::filter()  masks stats::filter()
## ✖ dplyr::lag()     masks stats::lag()
## ✖ purrr::map()     masks maps::map()
## ✖ dplyr::select()  masks raster::select()
## ℹ Use the conflicted package (<http://conflicted.r-lib.org/>) to force all conflicts to become errors
library(tigris)
## To enable caching of data, set `options(tigris_use_cache = TRUE)`
## in your R script or .Rprofile.
options(tigris_use_cache = TRUE)

Read in the GBIF data points

gbif = read.csv("../Formatted.Data/gbif.filtered.csv")

Subset for Abies balsamea

ABBA.occ = gbif %>%
  filter(species == "Abies balsamea") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

ABBA.range = st_read("../../USTreeAtlas/shp/abiebals/")
## Reading layer `abiebals' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/abiebals' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 442 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -116.6732 ymin: 38.40429 xmax: -52.61445 ymax: 58.79789
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(ABBA.range)  <- 4267 
maps::map(database = "state")
plot(ABBA.range$geometry, col="red", add=T)

us_states <- states(cb = TRUE)
## Retrieving data for the year 2021
continental_states <- us_states %>%
  filter(!NAME %in% (c("Alaska","American Samoa","Guam","Commonwealth of the Northern Mariana Islands","Hawaii","United States Virgin Islands",
                       "Puerto Rico")))
states.map = continental_states %>%
  st_as_sf %>%
  st_transform(st_crs(ABBA.range))
ggplot(states.map) +
  geom_sf()

ggplot(ABBA.range) +
  geom_sf()

ggplot()+
  geom_sf(data = states.map)+
  geom_sf(data = ABBA.range, col = "red")

ABBA_clipped = st_intersection(ABBA.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = ABBA_clipped, col = "red")+
  theme_classic()

ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = ABBA_clipped, col = "red", linewidth = 2)+
  geom_point(data = ABBA.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/ABBA.range.pdf", width = 12, height = 8)

Repeat basics of steps above for all other species

Subset for Abies fraseri

ABFR.occ = gbif %>%
  filter(species == "Abies fraseri") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

ABFR.range = st_read("../../USTreeAtlas/shp/abiefras/")
## Reading layer `abiefras' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/abiefras' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 8 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -83.82256 ymin: 35.03492 xmax: -81.02652 ymax: 37.00807
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(ABFR.range)  <- 4267 
ABFR_clipped = st_intersection(ABFR.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = ABFR_clipped, col = "red", linewidth = 2)+
  geom_point(data = ABFR.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/ABFR.range.pdf", width = 12, height = 8)

Subset for Acer barbatum

ACBA3.occ = gbif %>%
  filter(species == "Acer barbatum") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

ACBA3.range = st_read("../../USTreeAtlas/shp/acerbarb/")
## Reading layer `acerbarb' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/acerbarb' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 21 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -98.45653 ymin: 28.42883 xmax: -76.26695 ymax: 37.3846
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(ACBA3.range)  <- 4267 
ACBA3_clipped = st_intersection(ACBA3.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = ACBA3_clipped, col = "red", linewidth = 2)+
  geom_point(data = ACBA3.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/ACBA3.range.pdf", width = 12, height = 8)

Subset for Acer leucoderme

ACLE.occ = gbif %>%
  filter(species == "Acer leucoderme") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

ACLE.range = st_read("../../USTreeAtlas/shp/acerleuc/")
## Reading layer `acerleuc' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/acerleuc' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 21 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -95.10542 ymin: 30.2384 xmax: -78.57996 ymax: 36.43029
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(ACLE.range)  <- 4267 
ACLE_clipped = st_intersection(ACLE.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = ACLE_clipped, col = "red", linewidth = 2)+
  geom_point(data = ACLE.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/ACLE.range.pdf", width = 12, height = 8)

Subset for Acer negundo

ACNE3.occ = gbif %>%
  filter(species == "Acer negundo") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

ACNE3.range = st_read("../../USTreeAtlas/shp/acernegu/")
## Reading layer `acernegu' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/acernegu' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 114 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -123.3583 ymin: 28.30204 xmax: -67.68153 ymax: 54.28115
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(ACNE3.range)  <- 4267 
ACNE3_clipped = st_intersection(ACNE3.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = ACNE3_clipped, col = "red", linewidth = 2)+
  geom_point(data = ACNE3.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/ACNE3.range.pdf", width = 12, height = 8)

Subset for Acer nigrum

ACNI5.occ = gbif %>%
  filter(species == "Acer nigrum") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

ACNI5.range = st_read("../../USTreeAtlas/shp/acernigr/")
## Reading layer `acernigr' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/acernigr' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 92 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -96.94963 ymin: 35.15816 xmax: -72.13584 ymax: 45.95228
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(ACNI5.range)  <- 4267 
ACNI5_clipped = st_intersection(ACNI5.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = ACNI5_clipped, col = "red", linewidth = 2)+
  geom_point(data = ACNI5.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/ACNI5.range.pdf", width = 12, height = 8)

Subset for Acer pensylvanicum

ACPE.occ = gbif %>%
  filter(species == "Acer pensylvanicum") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

ACPE.range = st_read("../../USTreeAtlas/shp/acerpens/")
## Reading layer `acerpens' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/acerpens' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 139 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -89.90945 ymin: 34.813 xmax: -59.68972 ymax: 49.26055
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(ACPE.range)  <- 4267 
ACPE_clipped = st_intersection(ACPE.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = ACPE_clipped, col = "red", linewidth = 2)+
  geom_point(data = ACPE.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/ACPE.range.pdf", width = 12, height = 8)

Subset for Acer rubrum

ACRU.occ = gbif %>%
  filter(species == "Acer rubrum") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

ACRU.range = st_read("../../USTreeAtlas/shp/acerrubr/")
## Reading layer `acerrubr' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/acerrubr' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 498 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -97.25412 ymin: 25.5467 xmax: -52.61445 ymax: 49.401
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(ACRU.range)  <- 4267 
ACRU_clipped = st_intersection(ACRU.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = ACRU_clipped, col = "red", linewidth = 2)+
  geom_point(data = ACRU.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/ACRU.range.pdf", width = 12, height = 8)

Subset for Acer saccharinum

ACSA2.occ = gbif %>%
  filter(species == "Acer saccharinum") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

ACSA2.range = st_read("../../USTreeAtlas/shp/acersacc/")
## Reading layer `acersacc' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/acersacc' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 38 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -98.05198 ymin: 29.9023 xmax: -66.47555 ymax: 48.20345
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(ACSA2.range)  <- 4267 
ACSA2_clipped = st_intersection(ACSA2.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = ACSA2_clipped, col = "red", linewidth = 2)+
  geom_point(data = ACSA2.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/ACSA2.range.pdf", width = 12, height = 8)

Subset for Acer saccharum

ACSA3.occ = gbif %>%
  filter(species == "Acer saccharum") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

ACSA3.range = st_read("../../USTreeAtlas/shp/acersacr/")
## Reading layer `acersacr' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/acersacr' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 441 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -97.7326 ymin: 33.45683 xmax: -59.68972 ymax: 49.4221
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(ACSA3.range)  <- 4267 
ACSA3_clipped = st_intersection(ACSA3.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = ACSA3_clipped, col = "red", linewidth = 2)+
  geom_point(data = ACSA3.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/ACSA3.range.pdf", width = 12, height = 8)

Subset for Acer spicatum

ACSP2.occ = gbif %>%
  filter(species == "Acer spicatum") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

ACSP2.range = st_read("../../USTreeAtlas/shp/acerspic/")
## Reading layer `acerspic' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/acerspic' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 264 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -104.7974 ymin: 34.86549 xmax: -52.61445 ymax: 54.83275
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(ACSP2.range)  <- 4267 
ACSP2_clipped = st_intersection(ACSP2.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = ACSP2_clipped, col = "red", linewidth = 2)+
  geom_point(data = ACSP2.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/ACSP2.range.pdf", width = 12, height = 8)

Subset for Aesculus flava. Range map not available.

AEFL.occ = gbif %>%
  filter(species == "Aesculus flava") %>%
  select(species,decimalLatitude,decimalLongitude)
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_point(data = AEFL.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/AEFL.range.pdf", width = 12, height = 8)

Subset for Aesculus glabra

AEGL.occ = gbif %>%
  filter(species == "Aesculus glabra") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

AEGL.range = st_read("../../USTreeAtlas/shp/aescglab/")
## Reading layer `aescglab' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/aescglab' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 17 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -99.97602 ymin: 29.33802 xmax: -79.56816 ymax: 42.69205
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(AEGL.range)  <- 4267 
AEGL_clipped = st_intersection(AEGL.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = AEGL_clipped, col = "red", linewidth = 2)+
  geom_point(data = AEGL.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/AEGL.range.pdf", width = 12, height = 8)

Subset for Ailanthus altissima. Range map not available.

AIAL.occ = gbif %>%
  filter(species == "Ailanthus altissima") %>%
  select(species,decimalLatitude,decimalLongitude)
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_point(data = AIAL.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/AIAL.range.pdf", width = 12, height = 8)

Subset for Annona glabra. Range map not available.

ANGL4.occ = gbif %>%
  filter(species == "Annona glabra") %>%
  select(species,decimalLatitude,decimalLongitude)
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_point(data = ANGL4.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/ANGL4.range.pdf", width = 12, height = 8)

Subset for Asimina triloba

ASTR.occ = gbif %>%
  filter(species == "Asimina triloba") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

ASTR.range = st_read("../../USTreeAtlas/shp/asimtril/")
## Reading layer `asimtril' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/asimtril' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 44 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -97.22963 ymin: 29.65825 xmax: -74.09969 ymax: 43.43197
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(ASTR.range)  <- 4267 
ASTR_clipped = st_intersection(ASTR.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = ASTR_clipped, col = "red", linewidth = 2)+
  geom_point(data = ASTR.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/ASTR.range.pdf", width = 12, height = 8)

Subset for Betula alleghaniensis

BEAL2.occ = gbif %>%
  filter(species == "Betula alleghaniensis") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

BEAL2.range = st_read("../../USTreeAtlas/shp/betualle/")
## Reading layer `betualle' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/betualle' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 91 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -96.07653 ymin: 34.77789 xmax: -52.61445 ymax: 49.95027
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(BEAL2.range)  <- 4267 
BEAL2_clipped = st_intersection(BEAL2.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = BEAL2_clipped, col = "red", linewidth = 2)+
  geom_point(data = BEAL2.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/BEAL2.range.pdf", width = 12, height = 8)

Subset for Betula lenta

BELE.occ = gbif %>%
  filter(species == "Betula lenta") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

BELE.range = st_read("../../USTreeAtlas/shp/betulent/")
## Reading layer `betulent' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/betulent' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 20 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -87.99353 ymin: 33.1902 xmax: -68.21455 ymax: 48.70293
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(BELE.range)  <- 4267 
BELE_clipped = st_intersection(BELE.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = BELE_clipped, col = "red", linewidth = 2)+
  geom_point(data = BELE.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/BELE.range.pdf", width = 12, height = 8)

Subset for Betula nigra

BENI.occ = gbif %>%
  filter(species == "Betula nigra") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

BENI.range = st_read("../../USTreeAtlas/shp/betunigr/")
## Reading layer `betunigr' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/betunigr' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 48 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -96.42217 ymin: 29.56167 xmax: -70.57861 ymax: 45.09521
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(BENI.range)  <- 4267 
BENI_clipped = st_intersection(BENI.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = BENI_clipped, col = "red", linewidth = 2)+
  geom_point(data = BENI.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/BENI.range.pdf", width = 12, height = 8)

Subset for Betula papyrifera

BEPA.occ = gbif %>%
  filter(species == "Betula papyrifera") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

BEPA.range = st_read("../../USTreeAtlas/shp/betupapy/")
## Reading layer `betupapy' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/betupapy' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 482 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -164.7386 ymin: 35.027 xmax: -52.61445 ymax: 68.50498
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(BEPA.range)  <- 4267 
BEPA_clipped = st_intersection(BEPA.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = BEPA_clipped, col = "red", linewidth = 2)+
  geom_point(data = BEPA.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/BEPA.range.pdf", width = 12, height = 8)

Subset for Betula populifolia

BEPO.occ = gbif %>%
  filter(species == "Betula populifolia") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

BEPO.range = st_read("../../USTreeAtlas/shp/betupopu/")
## Reading layer `betupopu' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/betupopu' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 145 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -87.36114 ymin: 35.00822 xmax: -59.68972 ymax: 47.98228
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(BEPO.range)  <- 4267 
BEPO_clipped = st_intersection(BEPO.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = BEPO_clipped, col = "red", linewidth = 2)+
  geom_point(data = BEPO.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/BEPO.range.pdf", width = 12, height = 8)

Subset for Carpinus caroliniana

CACA18.occ = gbif %>%
  filter(species == "Carpinus caroliniana") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

CACA18.range = st_read("../../USTreeAtlas/shp/carpcaro/")
## Reading layer `carpcaro' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/carpcaro' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 74 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -108.0685 ymin: 13.74625 xmax: -68.47728 ymax: 47.47089
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(CACA18.range)  <- 4267 
CACA18_clipped = st_intersection(CACA18.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = CACA18_clipped, col = "red", linewidth = 2)+
  geom_point(data = CACA18.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/CACA18.range.pdf", width = 12, height = 8)

Subset for Carya alba. No range map available

CAAL27.occ = gbif %>%
  filter(species == "Carya alba") %>%
  select(species,decimalLatitude,decimalLongitude)
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_point(data = CAAL27.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/CAAL27.range.pdf", width = 12, height = 8)

Subset for Carya aquatica

CAAQ2.occ = gbif %>%
  filter(species == "Carya aquatica") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

CAAQ2.range = st_read("../../USTreeAtlas/shp/caryaqua/")
## Reading layer `caryaqua' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/caryaqua' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 46 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -97.70158 ymin: 26.41639 xmax: -75.71973 ymax: 37.88053
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(CAAQ2.range)  <- 4267 
CAAQ2_clipped = st_intersection(CAAQ2.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = CAAQ2_clipped, col = "red", linewidth = 2)+
  geom_point(data = CAAQ2.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/CAAQ2.range.pdf", width = 12, height = 8)

Subset for Carya cordiformis

CACO15.occ = gbif %>%
  filter(species == "Carya cordiformis") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

CACO15.range = st_read("../../USTreeAtlas/shp/carycord/")
## Reading layer `carycord' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/carycord' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 62 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -98.77852 ymin: 29.61891 xmax: -70.9063 ymax: 47.43512
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(CACO15.range)  <- 4267 
CACO15_clipped = st_intersection(CACO15.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = CACO15_clipped, col = "red", linewidth = 2)+
  geom_point(data = CACO15.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/CACO15.range.pdf", width = 12, height = 8)

Subset for Carya glabra

CAGL8.occ = gbif %>%
  filter(species == "Carya glabra") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

CAGL8.range = st_read("../../USTreeAtlas/shp/caryglab/")
## Reading layer `caryglab' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/caryglab' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 88 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -95.0925 ymin: 26.5667 xmax: -69.92722 ymax: 44.09084
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(CAGL8.range)  <- 4267 
CAGL8_clipped = st_intersection(CAGL8.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = CAGL8_clipped, col = "red", linewidth = 2)+
  geom_point(data = CAGL8.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/CAGL8.range.pdf", width = 12, height = 8)

Subset for Carya illinoinensis

CAIL2.occ = gbif %>%
  filter(species == "Carya illinoinensis") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

CAIL2.range = st_read("../../USTreeAtlas/shp/caryilli/")
## Reading layer `caryilli' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/caryilli' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 22 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -104.0874 ymin: 25.38815 xmax: -84.48216 ymax: 42.30187
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(CAIL2.range)  <- 4267 
CAIL2_clipped = st_intersection(CAIL2.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = CAIL2_clipped, col = "red", linewidth = 2)+
  geom_point(data = CAIL2.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/CAIL2.range.pdf", width = 12, height = 8)

Subset for Carya ovata

CAOV2.occ = gbif %>%
  filter(species == "Carya ovata") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

CAOV2.range = st_read("../../USTreeAtlas/shp/caryovat/")
## Reading layer `caryovat' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/caryovat' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 49 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -101.1183 ymin: 19.93666 xmax: -69.77362 ymax: 46.39629
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(CAOV2.range)  <- 4267 
CAOV2_clipped = st_intersection(CAOV2.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = CAOV2_clipped, col = "red", linewidth = 2)+
  geom_point(data = CAOV2.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/CAOV2.range.pdf", width = 12, height = 8)

Subset for Carya texana

CATE9.occ = gbif %>%
  filter(species == "Carya texana") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

CATE9.range = st_read("../../USTreeAtlas/shp/carytexa/")
## Reading layer `carytexa' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/carytexa' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 10 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -99.2131 ymin: 28.6377 xmax: -87.35191 ymax: 40.45395
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(CATE9.range)  <- 4267 
CATE9_clipped = st_intersection(CATE9.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = CATE9_clipped, col = "red", linewidth = 2)+
  geom_point(data = CATE9.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/CATE9.range.pdf", width = 12, height = 8)

Subset for Castanea dentata

CADE12.occ = gbif %>%
  filter(species == "Castanea dentata") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

CADE12.range = st_read("../../USTreeAtlas/shp/castdent/")
## Reading layer `castdent' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/castdent' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 69 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -90.06321 ymin: 30.78233 xmax: -68.80517 ymax: 44.94212
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(CADE12.range)  <- 4267 
CADE12_clipped = st_intersection(CADE12.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = CADE12_clipped, col = "red", linewidth = 2)+
  geom_point(data = CADE12.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/CADE12.range.pdf", width = 12, height = 8)

Subset for Celtis laevigata

CELA.occ = gbif %>%
  filter(species == "Celtis laevigata") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

CELA.range = st_read("../../USTreeAtlas/shp/celtlaev/")
## Reading layer `celtlaev' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/celtlaev' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 83 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -104.1844 ymin: 21.94096 xmax: -75.44749 ymax: 40.79288
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(CELA.range)  <- 4267 
CELA_clipped = st_intersection(CELA.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = CELA_clipped, col = "red", linewidth = 2)+
  geom_point(data = CELA.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/CELA.range.pdf", width = 12, height = 8)

Subset for Celtis occidentalis

CEOC.occ = gbif %>%
  filter(species == "Celtis occidentalis") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

CEOC.range = st_read("../../USTreeAtlas/shp/celtocci/")
## Reading layer `celtocci' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/celtocci' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 63 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -104.591 ymin: 32.21578 xmax: -71.17735 ymax: 50.1872
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(CEOC.range)  <- 4267 
CEOC_clipped = st_intersection(CEOC.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = CEOC_clipped, col = "red", linewidth = 2)+
  geom_point(data = CEOC.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/CEOC.range.pdf", width = 12, height = 8)

Subset for Cercis canadensis. Range map unavailable

CECA4.occ = gbif %>%
  filter(species == "Cercis canadensis") %>%
  select(species,decimalLatitude,decimalLongitude)
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_point(data = CECA4.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/CECA4.range.pdf", width = 12, height = 8)

Subset for Chamaecyparis thyoides

CHTH2.occ = gbif %>%
  filter(species == "Chamaecyparis thyoides") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

CHTH2.range = st_read("../../USTreeAtlas/shp/chamthyo/")
## Reading layer `chamthyo' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/chamthyo' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 67 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -89.87479 ymin: 29.12146 xmax: -68.76874 ymax: 44.42347
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(CHTH2.range)  <- 4267 
CHTH2_clipped = st_intersection(CHTH2.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = CHTH2_clipped, col = "red", linewidth = 2)+
  geom_point(data = CHTH2.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/CHTH2.range.pdf", width = 12, height = 8)

Subset for Cornus florida

COFL2.occ = gbif %>%
  filter(species == "Cornus florida") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

COFL2.range = st_read("../../USTreeAtlas/shp/cornflor/")
## Reading layer `cornflor' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/cornflor' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 84 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -100.8646 ymin: 18.52992 xmax: -69.92722 ymax: 43.94442
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(COFL2.range)  <- 4267 
COFL2_clipped = st_intersection(COFL2.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = COFL2_clipped, col = "red", linewidth = 2)+
  geom_point(data = COFL2.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/COFL2.range.pdf", width = 12, height = 8)

Subset for Diospyros virginiana

DIVI5.occ = gbif %>%
  filter(species == "Diospyros virginiana") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

DIVI5.range = st_read("../../USTreeAtlas/shp/diosvirg/")
## Reading layer `diosvirg' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/diosvirg' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 100 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -99.25706 ymin: 25.02361 xmax: -72.77579 ymax: 41.37863
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(DIVI5.range)  <- 4267 
DIVI5_clipped = st_intersection(DIVI5.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = DIVI5_clipped, col = "red", linewidth = 2)+
  geom_point(data = DIVI5.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/DIVI5.range.pdf", width = 12, height = 8)

Subset for Fagus grandifolia

FAGR.occ = gbif %>%
  filter(species == "Fagus grandifolia") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

FAGR.range = st_read("../../USTreeAtlas/shp/fagugran/")
## Reading layer `fagugran' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/fagugran' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 171 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -99.64175 ymin: 19.44994 xmax: -59.68972 ymax: 47.92627
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(FAGR.range)  <- 4267 
FAGR_clipped = st_intersection(FAGR.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = FAGR_clipped, col = "red", linewidth = 2)+
  geom_point(data = FAGR.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/FAGR.range.pdf", width = 12, height = 8)

Subset for Fraxinus americana

FRAM2.occ = gbif %>%
  filter(species == "Fraxinus americana") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

FRAM2.range = st_read("../../USTreeAtlas/shp/fraxamer/")
## Reading layer `fraxamer' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/fraxamer' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 157 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -98.82271 ymin: 28.62632 xmax: -59.79973 ymax: 48.1909
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(FRAM2.range)  <- 4267 
FRAM2_clipped = st_intersection(FRAM2.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = FRAM2_clipped, col = "red", linewidth = 2)+
  geom_point(data = FRAM2.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/FRAM2.range.pdf", width = 12, height = 8)

Subset for Fraxinus caroliniana

FRCA3.occ = gbif %>%
  filter(species == "Fraxinus caroliniana") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

FRCA3.range = st_read("../../USTreeAtlas/shp/fraxcaro/")
## Reading layer `fraxcaro' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/fraxcaro' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 71 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -95.62966 ymin: 25.74786 xmax: -75.44749 ymax: 38.66853
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(FRCA3.range)  <- 4267 
FRCA3_clipped = st_intersection(FRCA3.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = FRCA3_clipped, col = "red", linewidth = 2)+
  geom_point(data = FRCA3.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/FRCA3.range.pdf", width = 12, height = 8)

Subset for Fraxinus nigra

FRNI.occ = gbif %>%
  filter(species == "Fraxinus nigra") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

FRNI.range = st_read("../../USTreeAtlas/shp/fraxnigr/")
## Reading layer `fraxnigr' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/fraxnigr' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 144 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -97.86788 ymin: 37.5547 xmax: -57.8979 ymax: 51.50689
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(FRNI.range)  <- 4267 
FRNI_clipped = st_intersection(FRNI.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = FRNI_clipped, col = "red", linewidth = 2)+
  geom_point(data = FRNI.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/FRNI.range.pdf", width = 12, height = 8)

Subset for Fraxinus pennsylvanica

FRPE.occ = gbif %>%
  filter(species == "Fraxinus pennsylvanica") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

FRPE.range = st_read("../../USTreeAtlas/shp/fraxpenn/")
## Reading layer `fraxpenn' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/fraxpenn' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 188 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -111.6724 ymin: 28.31472 xmax: -59.68972 ymax: 54.27539
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(FRPE.range)  <- 4267 
FRPE_clipped = st_intersection(FRPE.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = FRPE_clipped, col = "red", linewidth = 2)+
  geom_point(data = FRPE.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/FRPE.range.pdf", width = 12, height = 8)

Subset for Fraxinus profunda

FRPR.occ = gbif %>%
  filter(species == "Fraxinus profunda") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

FRPR.range = st_read("../../USTreeAtlas/shp/fraxprof/")
## Reading layer `fraxprof' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/fraxprof' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 58 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -93.60954 ymin: 28.59934 xmax: -75.44749 ymax: 41.38225
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(FRPR.range)  <- 4267 
FRPR_clipped = st_intersection(FRPR.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = FRPR_clipped, col = "red", linewidth = 2)+
  geom_point(data = FRPR.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/FRPR.range.pdf", width = 12, height = 8)

Subset for Fraxinus quadrangulata

FRQU.occ = gbif %>%
  filter(species == "Fraxinus quadrangulata") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

FRQU.range = st_read("../../USTreeAtlas/shp/fraxquad/")
## Reading layer `fraxquad' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/fraxquad' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 18 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -96.2248 ymin: 33.2491 xmax: -81.49872 ymax: 42.96199
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(FRQU.range)  <- 4267 
FRQU_clipped = st_intersection(FRQU.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = FRQU_clipped, col = "red", linewidth = 2)+
  geom_point(data = FRQU.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/FRQU.range.pdf", width = 12, height = 8)

Subset for Fraxinus texensis/Fraxinus albicans (gbif)

FRTE.occ = gbif %>%
  filter(species == "Fraxinus albicans") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

FRTE.range = st_read("../../USTreeAtlas/shp/fraxtexe/")
## Reading layer `fraxtexe' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/fraxtexe' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 8 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -101.2223 ymin: 28.28765 xmax: -96.66921 ymax: 34.50977
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(FRTE.range)  <- 4267 
FRTE_clipped = st_intersection(FRTE.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = FRTE_clipped, col = "red", linewidth = 2)+
  geom_point(data = FRTE.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/FRTE.range.pdf", width = 12, height = 8)

Subset for Gleditsia aquatica

GLAQ.occ = gbif %>%
  filter(species == "Gleditsia aquatica") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

GLAQ.range = st_read("../../USTreeAtlas/shp/gledaqua/")
## Reading layer `gledaqua' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/gledaqua' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 20 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -96.36879 ymin: 27.27744 xmax: -79.58623 ymax: 38.97958
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(GLAQ.range)  <- 4267 
GLAQ_clipped = st_intersection(GLAQ.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = GLAQ_clipped, col = "red", linewidth = 2)+
  geom_point(data = GLAQ.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/GLAQ.range.pdf", width = 12, height = 8)

Subset for Gleditsia triacanthos

GLTR.occ = gbif %>%
  filter(species == "Gleditsia triacanthos") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

GLTR.range = st_read("../../USTreeAtlas/shp/gledtria/")
## Reading layer `gledtria' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/gledtria' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 18 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -101.5966 ymin: 28.73964 xmax: -76.56055 ymax: 43.67306
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(GLTR.range)  <- 4267 
GLTR_clipped = st_intersection(GLTR.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = GLTR_clipped, col = "red", linewidth = 2)+
  geom_point(data = GLTR.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/GLTR.range.pdf", width = 12, height = 8)

Subset for Gordonia lasianthus

GOLA.occ = gbif %>%
  filter(species == "Gordonia lasianthus") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

GOLA.range = st_read("../../USTreeAtlas/shp/gordlasi/")
## Reading layer `gordlasi' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/gordlasi' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 29 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -89.66769 ymin: 27.00125 xmax: -75.71973 ymax: 35.99722
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(GOLA.range)  <- 4267 
GOLA_clipped = st_intersection(GOLA.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = GOLA_clipped, col = "red", linewidth = 2)+
  geom_point(data = GOLA.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/GOLA.range.pdf", width = 12, height = 8)

Subset for Gymnocladus dioicus

GYDI.occ = gbif %>%
  filter(species == "Gymnocladus dioicus") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

GYDI.range = st_read("../../USTreeAtlas/shp/gymndioi/")
## Reading layer `gymndioi' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/gymndioi' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 44 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -99.97971 ymin: 33.98936 xmax: -76.18254 ymax: 45.07224
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(GYDI.range)  <- 4267 
GYDI_clipped = st_intersection(GYDI.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = GYDI_clipped, col = "red", linewidth = 2)+
  geom_point(data = GYDI.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/GYDI.range.pdf", width = 12, height = 8)

Subset for Ilex opaca

ILOP.occ = gbif %>%
  filter(species == "Ilex opaca") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

ILOP.range = st_read("../../USTreeAtlas/shp/ilexopac/")
## Reading layer `ilexopac' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/ilexopac' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 67 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -98.05645 ymin: 27.40893 xmax: -70.25293 ymax: 42.31931
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(ILOP.range)  <- 4267 
ILOP_clipped = st_intersection(ILOP.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = ILOP_clipped, col = "red", linewidth = 2)+
  geom_point(data = ILOP.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/ILOP.range.pdf", width = 12, height = 8)

Subset for Juglans cinerea

JUCI.occ = gbif %>%
  filter(species == "Juglans cinerea") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

JUCI.range = st_read("../../USTreeAtlas/shp/juglcine/")
## Reading layer `juglcine' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/juglcine' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 44 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -94.84976 ymin: 33.36165 xmax: -65.67808 ymax: 47.24971
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(JUCI.range)  <- 4267 
JUCI_clipped = st_intersection(JUCI.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = JUCI_clipped, col = "red", linewidth = 2)+
  geom_point(data = JUCI.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/JUCI.range.pdf", width = 12, height = 8)

Subset for Juglans nigra

JUNI.occ = gbif %>%
  filter(species == "Juglans nigra") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

JUNI.range = st_read("../../USTreeAtlas/shp/juglnigr/")
## Reading layer `juglnigr' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/juglnigr' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 61 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -100.5016 ymin: 28.65753 xmax: -71.5475 ymax: 45.06173
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(JUNI.range)  <- 4267 
JUNI_clipped = st_intersection(JUNI.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = JUNI_clipped, col = "red", linewidth = 2)+
  geom_point(data = JUNI.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/JUNI.range.pdf", width = 12, height = 8)

Subset for Juniperus virginiana

JUVI.occ = gbif %>%
  filter(species == "Juniperus virginiana") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

JUVI.range = st_read("../../USTreeAtlas/shp/junivirg/")
## Reading layer `junivirg' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/junivirg' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 60 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -103.5662 ymin: 29.45831 xmax: -69.77362 ymax: 47.07206
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(JUVI.range)  <- 4267 
JUVI_clipped = st_intersection(JUVI.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = JUVI_clipped, col = "red", linewidth = 2)+
  geom_point(data = JUVI.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/JUVI.range.pdf", width = 12, height = 8)

Subset for Larix laricina

LALA.occ = gbif %>%
  filter(species == "Larix laricina") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

LALA.range = st_read("../../USTreeAtlas/shp/larilari/")
## Reading layer `larilari' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/larilari' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 471 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -160.2373 ymin: 39.14292 xmax: -52.61445 ymax: 67.67176
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(LALA.range)  <- 4267 
LALA_clipped = st_intersection(LALA.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = LALA_clipped, col = "red", linewidth = 2)+
  geom_point(data = LALA.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/LALA.range.pdf", width = 12, height = 8)

Subset for Liquidambar styraciflua

LIST2.occ = gbif %>%
  filter(species == "Liquidambar styraciflua") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

LIST2.range = st_read("../../USTreeAtlas/shp/liqustyr/")
## Reading layer `liqustyr' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/liqustyr' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 81 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -99.09172 ymin: 13.20806 xmax: -73.26896 ymax: 41.15219
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(LIST2.range)  <- 4267 
LIST2_clipped = st_intersection(LIST2.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = LIST2_clipped, col = "red", linewidth = 2)+
  geom_point(data = LIST2.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/LIST2.range.pdf", width = 12, height = 8)

Subset for Liriodendron tulipifera

LITU.occ = gbif %>%
  filter(species == "Liriodendron tulipifera") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

LITU.range = st_read("../../USTreeAtlas/shp/lirituli/")
## Reading layer `lirituli' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/lirituli' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 61 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -93.14224 ymin: 28.39825 xmax: -71.36444 ymax: 44.53041
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(LITU.range)  <- 4267 
LITU_clipped = st_intersection(LITU.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = LITU_clipped, col = "red", linewidth = 2)+
  geom_point(data = LITU.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/LITU.range.pdf", width = 12, height = 8)

Subset for Maclura pomifera

MAPO.occ = gbif %>%
  filter(species == "Maclura pomifera") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

MAPO.range = st_read("../../USTreeAtlas/shp/maclpomi/")
## Reading layer `maclpomi' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/maclpomi' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 2 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -103.3863 ymin: 29.07236 xmax: -93.57389 ymax: 34.64141
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(MAPO.range)  <- 4267 
MAPO_clipped = st_intersection(MAPO.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = MAPO_clipped, col = "red", linewidth = 2)+
  geom_point(data = MAPO.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/MAPO.range.pdf", width = 12, height = 8)

Subset for Magnolia acuminata

MAAC.occ = gbif %>%
  filter(species == "Magnolia acuminata") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

MAAC.range = st_read("../../USTreeAtlas/shp/magnacum/")
## Reading layer `magnacum' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/magnacum' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 49 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -94.6578 ymin: 30.68451 xmax: -75.45156 ymax: 43.38474
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(MAAC.range)  <- 4267 
MAAC_clipped = st_intersection(MAAC.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = MAAC_clipped, col = "red", linewidth = 2)+
  geom_point(data = MAAC.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/MAAC.range.pdf", width = 12, height = 8)

Subset for Magnolia fraseri

MAFR.occ = gbif %>%
  filter(species == "Magnolia fraseri") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

MAFR.range = st_read("../../USTreeAtlas/shp/magnfras/")
## Reading layer `magnfras' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/magnfras' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 4 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -84.82449 ymin: 34.46459 xmax: -77.81695 ymax: 39.62026
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(MAFR.range)  <- 4267 
MAFR_clipped = st_intersection(MAFR.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = MAFR_clipped, col = "red", linewidth = 2)+
  geom_point(data = MAFR.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/MAFR.range.pdf", width = 12, height = 8)

Subset for Magnolia grandiflora

MAGR4.occ = gbif %>%
  filter(species == "Magnolia grandiflora") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

MAGR4.range = st_read("../../USTreeAtlas/shp/magngran/")
## Reading layer `magngran' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/magngran' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 28 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -96.0386 ymin: 27.16424 xmax: -76.68305 ymax: 35.66163
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(MAGR4.range)  <- 4267 
MAGR4_clipped = st_intersection(MAGR4.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = MAGR4_clipped, col = "red", linewidth = 2)+
  geom_point(data = MAGR4.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/MAGR4.range.pdf", width = 12, height = 8)

Subset for Magnolia macrophylla

MAMA2.occ = gbif %>%
  filter(species == "Magnolia macrophylla") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

MAMA2.range = st_read("../../USTreeAtlas/shp/magnmacr/")
## Reading layer `magnmacr' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/magnmacr' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 25 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -92.80724 ymin: 30.45176 xmax: -78.51603 ymax: 39.02176
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(MAMA2.range)  <- 4267 
MAMA2_clipped = st_intersection(MAMA2.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = MAMA2_clipped, col = "red", linewidth = 2)+
  geom_point(data = MAMA2.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/MAMA2.range.pdf", width = 12, height = 8)

Subset for Magnolia virginiana

MAVI2.occ = gbif %>%
  filter(species == "Magnolia virginiana") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

MAVI2.range = st_read("../../USTreeAtlas/shp/magnvirg/")
## Reading layer `magnvirg' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/magnvirg' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 82 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -95.6114 ymin: 25.71463 xmax: -72.75641 ymax: 40.91171
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(MAVI2.range)  <- 4267 
MAVI2_clipped = st_intersection(MAVI2.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = MAVI2_clipped, col = "red", linewidth = 2)+
  geom_point(data = MAVI2.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/MAVI2.range.pdf", width = 12, height = 8)

Subset for Morus rubra

MORU2.occ = gbif %>%
  filter(species == "Morus rubra") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

MORU2.range = st_read("../../USTreeAtlas/shp/morurubr/")
## Reading layer `morurubr' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/morurubr' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 96 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -101.9412 ymin: 25.71206 xmax: -70.92392 ymax: 44.91562
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(MORU2.range)  <- 4267 
MORU2_clipped = st_intersection(MORU2.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = MORU2_clipped, col = "red", linewidth = 2)+
  geom_point(data = MORU2.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/MORU2.range.pdf", width = 12, height = 8)

Subset for Nyssa aquatica

NYAQ2.occ = gbif %>%
  filter(species == "Nyssa aquatica") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

NYAQ2.range = st_read("../../USTreeAtlas/shp/nyssaqua/")
## Reading layer `nyssaqua' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/nyssaqua' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 48 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -95.18513 ymin: 29.33239 xmax: -75.44749 ymax: 39.09947
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(NYAQ2.range)  <- 4267 
NYAQ2_clipped = st_intersection(NYAQ2.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = NYAQ2_clipped, col = "red", linewidth = 2)+
  geom_point(data = NYAQ2.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/NYAQ2.range.pdf", width = 12, height = 8)

Subset for Nyssa biflora. No range map available.

NYBI.occ = gbif %>%
  filter(species == "Nyssa biflora") %>%
  select(species,decimalLatitude,decimalLongitude)
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_point(data = NYBI.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/NYBI.range.pdf", width = 12, height = 8)

Subset for Nyssa ogeche

NYOG.occ = gbif %>%
  filter(species == "Nyssa ogeche") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

NYOG.range = st_read("../../USTreeAtlas/shp/nyssogec/")
## Reading layer `nyssogec' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/nyssogec' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 8 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -86.51112 ymin: 29.52137 xmax: -80.47782 ymax: 32.7832
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(NYOG.range)  <- 4267 
NYOG_clipped = st_intersection(NYOG.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = NYOG_clipped, col = "red", linewidth = 2)+
  geom_point(data = NYOG.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/NYOG.range.pdf", width = 12, height = 8)

Subset for Nyssa sylvatica

NYSY.occ = gbif %>%
  filter(species == "Nyssa sylvatica") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

NYSY.range = st_read("../../USTreeAtlas/shp/nysssylv/")
## Reading layer `nysssylv' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/nysssylv' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 98 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -98.92374 ymin: 16.64896 xmax: -69.58422 ymax: 44.80475
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(NYSY.range)  <- 4267 
NYSY_clipped = st_intersection(NYSY.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = NYSY_clipped, col = "red", linewidth = 2)+
  geom_point(data = NYSY.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/NYSY.range.pdf", width = 12, height = 8)

Subset for Ostrya virginiana

OSVI.occ = gbif %>%
  filter(species == "Ostrya virginiana") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

OSVI.range = st_read("../../USTreeAtlas/shp/ostrvirg/")
## Reading layer `ostrvirg' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/ostrvirg' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 145 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -109.4911 ymin: 13.70159 xmax: -59.68972 ymax: 50.35658
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(OSVI.range)  <- 4267 
OSVI_clipped = st_intersection(OSVI.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = OSVI_clipped, col = "red", linewidth = 2)+
  geom_point(data = OSVI.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/OSVI.range.pdf", width = 12, height = 8)

Subset for Oxydendrum arboreum

OXAR.occ = gbif %>%
  filter(species == "Oxydendrum arboreum") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

OXAR.range = st_read("../../USTreeAtlas/shp/oxydarbo/")
## Reading layer `oxydarbo' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/oxydarbo' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 34 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -94.92362 ymin: 29.90913 xmax: -75.31862 ymax: 40.52448
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(OXAR.range)  <- 4267 
OXAR_clipped = st_intersection(OXAR.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = OXAR_clipped, col = "red", linewidth = 2)+
  geom_point(data = OXAR.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/OXAR.range.pdf", width = 12, height = 8)

Subset for Paulownia tomentosa. No range map available.

PATO2.occ = gbif %>%
  filter(species == "Paulownia tomentosa") %>%
  select(species,decimalLatitude,decimalLongitude)
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_point(data = PATO2.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/PATO2.range.pdf", width = 12, height = 8)

Subset for Persea borbonia

PEBO.occ = gbif %>%
  filter(species == "Persea borbonia") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

PEBO.range = st_read("../../USTreeAtlas/shp/persborb/")
## Reading layer `persborb' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/persborb' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 94 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -97.96335 ymin: 24.6522 xmax: -75.1954 ymax: 38.55726
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(PEBO.range)  <- 4267 
PEBO_clipped = st_intersection(PEBO.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = PEBO_clipped, col = "red", linewidth = 2)+
  geom_point(data = PEBO.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/PEBO.range.pdf", width = 12, height = 8)

Subset for Picea glauca

PIGL.occ = gbif %>%
  filter(species == "Picea glauca") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

PIGL.range = st_read("../../USTreeAtlas/shp/piceglau/")
## Reading layer `piceglau' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/piceglau' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 1449 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -163.9861 ymin: 42.71046 xmax: -52.63629 ymax: 69.57938
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(PIGL.range)  <- 4267 
PIGL_clipped = st_intersection(PIGL.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = PIGL_clipped, col = "red", linewidth = 2)+
  geom_point(data = PIGL.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/PIGL.range.pdf", width = 12, height = 8)

Subset for Picea mariana

PIMA.occ = gbif %>%
  filter(species == "Picea mariana") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

PIMA.range = st_read("../../USTreeAtlas/shp/picemari/")
## Reading layer `picemari' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/picemari' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 1585 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -161.9129 ymin: 40.52763 xmax: -52.63629 ymax: 69.51231
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(PIMA.range)  <- 4267 
PIMA_clipped = st_intersection(PIMA.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = PIMA_clipped, col = "red", linewidth = 2)+
  geom_point(data = PIMA.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/PIMA.range.pdf", width = 12, height = 8)

Subset for Picea rubens

PIRU.occ = gbif %>%
  filter(species == "Picea rubens") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

PIRU.range = st_read("../../USTreeAtlas/shp/picerube/")
## Reading layer `picerube' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/picerube' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 326 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -83.7767 ymin: 35.23503 xmax: -59.69256 ymax: 48.01405
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(PIRU.range)  <- 4267 
PIRU_clipped = st_intersection(PIRU.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = PIRU_clipped, col = "red", linewidth = 2)+
  geom_point(data = PIRU.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/PIRU.range.pdf", width = 12, height = 8)

Subset for Pinus banksiana

PIBA2.occ = gbif %>%
  filter(species == "Pinus banksiana") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

PIBA2.range = st_read("../../USTreeAtlas/shp/pinubank/")
## Reading layer `pinubank' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/pinubank' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 395 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -127.2683 ymin: 41.47121 xmax: -60.3597 ymax: 65.29585
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(PIBA2.range)  <- 4267 
PIBA2_clipped = st_intersection(PIBA2.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = PIBA2_clipped, col = "red", linewidth = 2)+
  geom_point(data = PIBA2.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/PIBA2.range.pdf", width = 12, height = 8)

Subset for Pinus clausa

PICL.occ = gbif %>%
  filter(species == "Pinus clausa") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

PICL.range = st_read("../../USTreeAtlas/shp/pinuclau/")
## Reading layer `pinuclau' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/pinuclau' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 52 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -87.70842 ymin: 25.84471 xmax: -80.03255 ymax: 30.66881
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(PICL.range)  <- 4267 
PICL_clipped = st_intersection(PICL.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = PICL_clipped, col = "red", linewidth = 2)+
  geom_point(data = PICL.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/PICL.range.pdf", width = 12, height = 8)

Subset for Pinus echinata

PIEC2.occ = gbif %>%
  filter(species == "Pinus echinata") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

PIEC2.range = st_read("../../USTreeAtlas/shp/pinuechi/")
## Reading layer `pinuechi' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/pinuechi' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 68 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -96.62461 ymin: 29.79749 xmax: -73.2604 ymax: 41.11497
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(PIEC2.range)  <- 4267 
PIEC2_clipped = st_intersection(PIEC2.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = PIEC2_clipped, col = "red", linewidth = 2)+
  geom_point(data = PIEC2.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/PIEC2.range.pdf", width = 12, height = 8)

Subset for Pinus elliottii

PIEL.occ = gbif %>%
  filter(species == "Pinus elliottii") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

PIEL.range = st_read("../../USTreeAtlas/shp/pinuelli/")
## Reading layer `pinuelli' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/pinuelli' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 129 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -90.77222 ymin: 24.58278 xmax: -79.60842 ymax: 33.28711
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(PIEL.range)  <- 4267 
PIEL_clipped = st_intersection(PIEL.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = PIEL_clipped, col = "red", linewidth = 2)+
  geom_point(data = PIEL.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/PIEL.range.pdf", width = 12, height = 8)

Subset for Pinus glabra

PIGL2.occ = gbif %>%
  filter(species == "Pinus glabra") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

PIGL2.range = st_read("../../USTreeAtlas/shp/pinuglab/")
## Reading layer `pinuglab' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/pinuglab' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 9 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -91.29132 ymin: 29.58212 xmax: -79.20943 ymax: 34.05908
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(PIGL2.range)  <- 4267 
PIGL2_clipped = st_intersection(PIGL2.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = PIGL2_clipped, col = "red", linewidth = 2)+
  geom_point(data = PIGL2.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/PIGL2.range.pdf", width = 12, height = 8)

Subset for Pinus palustris

PIPA2.occ = gbif %>%
  filter(species == "Pinus palustris") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

PIPA2.range = st_read("../../USTreeAtlas/shp/pinupalu/")
## Reading layer `pinupalu' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/pinupalu' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 38 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -95.21806 ymin: 26.61921 xmax: -75.79999 ymax: 36.84856
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(PIPA2.range)  <- 4267 
PIPA2_clipped = st_intersection(PIPA2.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = PIPA2_clipped, col = "red", linewidth = 2)+
  geom_point(data = PIPA2.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/PIPA2.range.pdf", width = 12, height = 8)

Subset for Pinus pungens

PIPUS.occ = gbif %>%
  filter(species == "Pinus pungens") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

PIPUS.range = st_read("../../USTreeAtlas/shp/pinupung/")
## Reading layer `pinupung' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/pinupung' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 26 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -84.29883 ymin: 34.56113 xmax: -74.77565 ymax: 41.06224
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(PIPUS.range)  <- 4267 
PIPUS_clipped = st_intersection(PIPUS.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = PIPUS_clipped, col = "red", linewidth = 2)+
  geom_point(data = PIPUS.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/PIPUS.range.pdf", width = 12, height = 8)

Subset for Pinus resinosa

PIRE.occ = gbif %>%
  filter(species == "Pinus resinosa") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

PIRE.range = st_read("../../USTreeAtlas/shp/pinuresi/")
## Reading layer `pinuresi' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/pinuresi' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 485 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -96.63834 ymin: 38.64299 xmax: -54.26736 ymax: 51.25153
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(PIRE.range)  <- 4267 
PIRE_clipped = st_intersection(PIRE.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = PIRE_clipped, col = "red", linewidth = 2)+
  geom_point(data = PIRE.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/PIRE.range.pdf", width = 12, height = 8)

Subset for Pinus rigida

PIRI.occ = gbif %>%
  filter(species == "Pinus rigida") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

PIRI.range = st_read("../../USTreeAtlas/shp/pinurigi/")
## Reading layer `pinurigi' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/pinurigi' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 154 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -88.39619 ymin: 33.68433 xmax: -68.12421 ymax: 45.25713
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(PIRI.range)  <- 4267 
PIRI_clipped = st_intersection(PIRI.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = PIRI_clipped, col = "red", linewidth = 2)+
  geom_point(data = PIRI.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/PIRI.range.pdf", width = 12, height = 8)

Subset for Pinus serotina

PISE.occ = gbif %>%
  filter(species == "Pinus serotina") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

PISE.range = st_read("../../USTreeAtlas/shp/pinusero/")
## Reading layer `pinusero' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/pinusero' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 40 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -87.52861 ymin: 28.02775 xmax: -74.69576 ymax: 39.87193
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(PISE.range)  <- 4267 
PISE_clipped = st_intersection(PISE.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = PISE_clipped, col = "red", linewidth = 2)+
  geom_point(data = PISE.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/PISE.range.pdf", width = 12, height = 8)

Subset for Pinus strobus

PIST.occ = gbif %>%
  filter(species == "Pinus strobus") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

PIST.range = st_read("../../USTreeAtlas/shp/pinustrb/")
## Reading layer `pinustrb' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/pinustrb' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 264 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -95.83807 ymin: 34.27946 xmax: -52.61445 ymax: 51.11618
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(PIST.range)  <- 4267 
PIST_clipped = st_intersection(PIST.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = PIST_clipped, col = "red", linewidth = 2)+
  geom_point(data = PIST.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/PIST.range.pdf", width = 12, height = 8)

Subset for Pinus taeda

PITA.occ = gbif %>%
  filter(species == "Pinus taeda") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

PITA.range = st_read("../../USTreeAtlas/shp/pinutaed/")
## Reading layer `pinutaed' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/pinutaed' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 58 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -97.6138 ymin: 28.22558 xmax: -75.01347 ymax: 39.50686
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(PITA.range)  <- 4267 
PITA_clipped = st_intersection(PITA.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = PITA_clipped, col = "red", linewidth = 2)+
  geom_point(data = PITA.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/PITA.range.pdf", width = 12, height = 8)

Subset for Pinus virginiana

PIVI2.occ = gbif %>%
  filter(species == "Pinus virginiana") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

PIVI2.range = st_read("../../USTreeAtlas/shp/pinuvirg/")
## Reading layer `pinuvirg' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/pinuvirg' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 36 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -89.41452 ymin: 32.30019 xmax: -73.15847 ymax: 41.37672
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(PIVI2.range)  <- 4267 
PIVI2_clipped = st_intersection(PIVI2.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = PIVI2_clipped, col = "red", linewidth = 2)+
  geom_point(data = PIVI2.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/PIVI2.range.pdf", width = 12, height = 8)

Subset for Planera aquatica

PLAQ.occ = gbif %>%
  filter(species == "Planera aquatica") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

PLAQ.range = st_read("../../USTreeAtlas/shp/planaqua/")
## Reading layer `planaqua' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/planaqua' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 17 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -96.84963 ymin: 29.49777 xmax: -78.27449 ymax: 37.5579
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(PLAQ.range)  <- 4267 
PLAQ_clipped = st_intersection(PLAQ.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = PLAQ_clipped, col = "red", linewidth = 2)+
  geom_point(data = PLAQ.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/PLAQ.range.pdf", width = 12, height = 8)

Subset for Platanus occidentalis

PLOC.occ = gbif %>%
  filter(species == "Platanus occidentalis") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

PLOC.range = st_read("../../USTreeAtlas/shp/platocci/")
## Reading layer `platocci' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/platocci' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 77 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -102.2446 ymin: 21.80153 xmax: -69.63906 ymax: 44.61452
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(PLOC.range)  <- 4267 
PLOC_clipped = st_intersection(PLOC.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = PLOC_clipped, col = "red", linewidth = 2)+
  geom_point(data = PLOC.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/PLOC.range.pdf", width = 12, height = 8)

Subset for Populus balsamifera

POBA2.occ = gbif %>%
  filter(species == "Populus balsamifera") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

POBA2.range = st_read("../../USTreeAtlas/shp/popubals/")
## Reading layer `popubals' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/popubals' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 415 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -163.9213 ymin: 38.58866 xmax: -52.61445 ymax: 69.24193
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(POBA2.range)  <- 4267 
POBA2_clipped = st_intersection(POBA2.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = POBA2_clipped, col = "red", linewidth = 2)+
  geom_point(data = POBA2.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/POBA2.range.pdf", width = 12, height = 8)

Subset for Populus deltoides

PODE3.occ = gbif %>%
  filter(species == "Populus deltoides") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

PODE3.range = st_read("../../USTreeAtlas/shp/popudelt/")
## Reading layer `popudelt' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/popudelt' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 79 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -114.3961 ymin: 27.32564 xmax: -72.24688 ymax: 52.14172
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(PODE3.range)  <- 4267 
PODE3_clipped = st_intersection(PODE3.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = PODE3_clipped, col = "red", linewidth = 2)+
  geom_point(data = PODE3.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/PODE3.range.pdf", width = 12, height = 8)

Subset for Populus grandidentata

POGR4.occ = gbif %>%
  filter(species == "Populus grandidentata") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

POGR4.range = st_read("../../USTreeAtlas/shp/popugran/")
## Reading layer `popugran' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/popugran' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 157 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -96.50261 ymin: 35.13372 xmax: -59.68972 ymax: 49.99322
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(POGR4.range)  <- 4267 
POGR4_clipped = st_intersection(POGR4.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = POGR4_clipped, col = "red", linewidth = 2)+
  geom_point(data = POGR4.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/POGR4.range.pdf", width = 12, height = 8)

Subset for Populus heterophylla

POHE4.occ = gbif %>%
  filter(species == "Populus heterophylla") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

POHE4.range = st_read("../../USTreeAtlas/shp/popuhete/")
## Reading layer `popuhete' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/popuhete' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 79 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -92.28237 ymin: 29.70388 xmax: -71.85611 ymax: 42.36729
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(POHE4.range)  <- 4267 
POHE4_clipped = st_intersection(POHE4.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = POHE4_clipped, col = "red", linewidth = 2)+
  geom_point(data = POHE4.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/POHE4.range.pdf", width = 12, height = 8)

Subset for Populus tremuloides

POTR5.occ = gbif %>%
  filter(species == "Populus tremuloides") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

POTR5.range = st_read("../../USTreeAtlas/shp/poputrem/")
## Reading layer `poputrem' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/poputrem' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 534 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -163.7564 ymin: 20.77411 xmax: -52.61445 ymax: 69.2738
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(POTR5.range)  <- 4267 
POTR5_clipped = st_intersection(POTR5.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = POTR5_clipped, col = "red", linewidth = 2)+
  geom_point(data = POTR5.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/POTR5.range.pdf", width = 12, height = 8)

Subset for Prunus americana

PRAM.occ = gbif %>%
  filter(species == "Prunus americana") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

PRAM.range = st_read("../../USTreeAtlas/shp/prunamer/")
## Reading layer `prunamer' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/prunamer' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 71 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -114.464 ymin: 28.57408 xmax: -70.17076 ymax: 50.03462
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(PRAM.range)  <- 4267 
PRAM_clipped = st_intersection(PRAM.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = PRAM_clipped, col = "red", linewidth = 2)+
  geom_point(data = PRAM.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/PRAM.range.pdf", width = 12, height = 8)

Subset for Prunus pensylvanica

PRPE2.occ = gbif %>%
  filter(species == "Prunus pensylvanica") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

PRPE2.range = st_read("../../USTreeAtlas/shp/prunpens/")
## Reading layer `prunpens' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/prunpens' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 336 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -128.2946 ymin: 34.80802 xmax: -52.61445 ymax: 62.31179
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(PRPE2.range)  <- 4267 
PRPE2_clipped = st_intersection(PRPE2.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = PRPE2_clipped, col = "red", linewidth = 2)+
  geom_point(data = PRPE2.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/PRPE2.range.pdf", width = 12, height = 8)

Subset for Prunus serotina

PRSE2.occ = gbif %>%
  filter(species == "Prunus serotina") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

PRSE2.range = st_read("../../USTreeAtlas/shp/prunsero/")
## Reading layer `prunsero' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/prunsero' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 200 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -113.9141 ymin: 14.56384 xmax: -63.3591 ymax: 47.8594
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(PRSE2.range)  <- 4267 
PRSE2_clipped = st_intersection(PRSE2.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = PRSE2_clipped, col = "red", linewidth = 2)+
  geom_point(data = PRSE2.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/PRSE2.range.pdf", width = 12, height = 8)

Subset for Prunus virginiana

PRVI.occ = gbif %>%
  filter(species == "Prunus virginiana") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

PRVI.range = st_read("../../USTreeAtlas/shp/prunvirg/")
## Reading layer `prunvirg' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/prunvirg' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 397 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -130.5808 ymin: 30.24095 xmax: -53.03568 ymax: 61.54309
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(PRVI.range)  <- 4267 
PRVI_clipped = st_intersection(PRVI.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = PRVI_clipped, col = "red", linewidth = 2)+
  geom_point(data = PRVI.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/PRVI.range.pdf", width = 12, height = 8)

Subset for Quercus alba

QUAL.occ = gbif %>%
  filter(species == "Quercus alba") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

QUAL.range = st_read("../../USTreeAtlas/shp/queralba/")
## Reading layer `queralba' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/queralba' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 82 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -96.34524 ymin: 29.63241 xmax: -69.11673 ymax: 46.45134
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(QUAL.range)  <- 4267 
QUAL_clipped = st_intersection(QUAL.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = QUAL_clipped, col = "red", linewidth = 2)+
  geom_point(data = QUAL.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/QUAL.range.pdf", width = 12, height = 8)

Subset for Quercus bicolor

QUBI.occ = gbif %>%
  filter(species == "Quercus bicolor") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

QUBI.range = st_read("../../USTreeAtlas/shp/querbico/")
## Reading layer `querbico' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/querbico' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 49 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -96.39155 ymin: 35.17383 xmax: -70.02106 ymax: 46.36058
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(QUBI.range)  <- 4267 
QUBI_clipped = st_intersection(QUBI.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = QUBI_clipped, col = "red", linewidth = 2)+
  geom_point(data = QUBI.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/QUBI.range.pdf", width = 12, height = 8)

Subset for Quercus coccinea

QUCO2.occ = gbif %>%
  filter(species == "Quercus coccinea") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

QUCO2.range = st_read("../../USTreeAtlas/shp/quercocc/")
## Reading layer `quercocc' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/quercocc' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 76 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -92.35469 ymin: 31.22379 xmax: -69.87 ymax: 44.1226
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(QUCO2.range)  <- 4267 
QUCO2_clipped = st_intersection(QUCO2.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = QUCO2_clipped, col = "red", linewidth = 2)+
  geom_point(data = QUCO2.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/QUCO2.range.pdf", width = 12, height = 8)

Subset for Quercus ellipsoidalis

QUEL.occ = gbif %>%
  filter(species == "Quercus ellipsoidalis") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

QUEL.range = st_read("../../USTreeAtlas/shp/querelli/")
## Reading layer `querelli' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/querelli' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 22 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -96.91765 ymin: 40.28243 xmax: -82.88901 ymax: 48.66541
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(QUEL.range)  <- 4267 
QUEL_clipped = st_intersection(QUEL.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = QUEL_clipped, col = "red", linewidth = 2)+
  geom_point(data = QUEL.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/QUEL.range.pdf", width = 12, height = 8)

Subset for Quercus falcata

QUFA.occ = gbif %>%
  filter(species == "Quercus falcata") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

QUFA.range = st_read("../../USTreeAtlas/shp/querfalc/")
## Reading layer `querfalc' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/querfalc' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 58 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -96.82243 ymin: 28.8895 xmax: -73.46515 ymax: 40.73988
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(QUFA.range)  <- 4267 
QUFA_clipped = st_intersection(QUFA.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = QUFA_clipped, col = "red", linewidth = 2)+
  geom_point(data = QUFA.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/QUFA.range.pdf", width = 12, height = 8)

Subset for Quercus ilicifolia

QUIL.occ = gbif %>%
  filter(species == "Quercus ilicifolia") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

QUIL.range = st_read("../../USTreeAtlas/shp/querilic/")
## Reading layer `querilic' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/querilic' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 46 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -81.30976 ymin: 35.22381 xmax: -67.91949 ymax: 44.84274
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(QUIL.range)  <- 4267 
QUIL_clipped = st_intersection(QUIL.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = QUIL_clipped, col = "red", linewidth = 2)+
  geom_point(data = QUIL.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/QUIL.range.pdf", width = 12, height = 8)

Subset for Quercus imbricaria

QUIM.occ = gbif %>%
  filter(species == "Quercus imbricaria") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

QUIM.range = st_read("../../USTreeAtlas/shp/querimbr/")
## Reading layer `querimbr' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/querimbr' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 38 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -96.96888 ymin: 31.03043 xmax: -75.56732 ymax: 42.31751
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(QUIM.range)  <- 4267 
QUIM_clipped = st_intersection(QUIM.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = QUIM_clipped, col = "red", linewidth = 2)+
  geom_point(data = QUIM.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/QUIM.range.pdf", width = 12, height = 8)

Subset for Quercus incana

QUIN.occ = gbif %>%
  filter(species == "Quercus incana") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

QUIN.range = st_read("../../USTreeAtlas/shp/querinca/")
## Reading layer `querinca' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/querinca' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 55 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -97.78163 ymin: 26.41639 xmax: -75.77622 ymax: 36.8999
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(QUIN.range)  <- 4267 
QUIN_clipped = st_intersection(QUIN.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = QUIN_clipped, col = "red", linewidth = 2)+
  geom_point(data = QUIN.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/QUIN.range.pdf", width = 12, height = 8)

Subset for Quercus laevis

QULA2.occ = gbif %>%
  filter(species == "Quercus laevis") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

QULA2.range = st_read("../../USTreeAtlas/shp/querlaev/")
## Reading layer `querlaev' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/querlaev' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 44 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -90.19388 ymin: 26.39276 xmax: -75.75282 ymax: 36.97728
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(QULA2.range)  <- 4267 
QULA2_clipped = st_intersection(QULA2.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = QULA2_clipped, col = "red", linewidth = 2)+
  geom_point(data = QULA2.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/QULA2.range.pdf", width = 12, height = 8)

Subset for Quercus laurifolia

QULA3.occ = gbif %>%
  filter(species == "Quercus laurifolia") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

QULA3.range = st_read("../../USTreeAtlas/shp/querlaur/")
## Reading layer `querlaur' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/querlaur' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 73 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -97.66097 ymin: 25.38946 xmax: -75.44749 ymax: 37.41531
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(QULA3.range)  <- 4267 
QULA3_clipped = st_intersection(QULA3.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = QULA3_clipped, col = "red", linewidth = 2)+
  geom_point(data = QULA3.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/QULA3.range.pdf", width = 12, height = 8)

Subset for Quercus lyrata

QULY.occ = gbif %>%
  filter(species == "Quercus lyrata") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

QULY.range = st_read("../../USTreeAtlas/shp/querlyra/")
## Reading layer `querlyra' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/querlyra' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 41 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -97.41548 ymin: 28.61861 xmax: -75.03833 ymax: 41.84388
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(QULY.range)  <- 4267 
QULY_clipped = st_intersection(QULY.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = QULY_clipped, col = "red", linewidth = 2)+
  geom_point(data = QULY.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/QULY.range.pdf", width = 12, height = 8)

Subset for Quercus macrocarpa

QUMA2.occ = gbif %>%
  filter(species == "Quercus macrocarpa") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

QUMA2.range = st_read("../../USTreeAtlas/shp/quermacr/")
## Reading layer `quermacr' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/quermacr' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 41 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -104.4497 ymin: 28.01439 xmax: -66.08747 ymax: 52.74889
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(QUMA2.range)  <- 4267 
QUMA2_clipped = st_intersection(QUMA2.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = QUMA2_clipped, col = "red", linewidth = 2)+
  geom_point(data = QUMA2.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/QUMA2.range.pdf", width = 12, height = 8)

Subset for Quercus margaretta. No range map available

QUMA6.occ = gbif %>%
  filter(species == "Quercus margaretta") %>%
  select(species,decimalLatitude,decimalLongitude)
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_point(data = QUMA6.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/QUMA6.range.pdf", width = 12, height = 8)

Subset for Quercus marilandica

QUMA3.occ = gbif %>%
  filter(species == "Quercus marilandica") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

QUMA3.range = st_read("../../USTreeAtlas/shp/quermari/")
## Reading layer `quermari' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/quermari' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 53 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -99.95763 ymin: 28.12027 xmax: -73.53086 ymax: 42.7558
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(QUMA3.range)  <- 4267 
QUMA3_clipped = st_intersection(QUMA3.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = QUMA3_clipped, col = "red", linewidth = 2)+
  geom_point(data = QUMA3.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/QUMA3.range.pdf", width = 12, height = 8)

Subset for Quercus michauxii (gbif) and Quercus prinus

QUMI.occ = gbif %>%
  filter(species == "Quercus michauxii") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

QUMI.range = st_read("../../USTreeAtlas/shp/quermich/")
## Reading layer `quermich' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/quermich' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 64 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -95.4984 ymin: 28.78528 xmax: -74.28231 ymax: 40.95467
## CRS:           NA
QUPR2.range = st_read("../../USTreeAtlas/shp/querprin/")
## Reading layer `querprin' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/querprin' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 29 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -89.98369 ymin: 31.95997 xmax: -70.50323 ymax: 44.57826
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(QUMI.range)  <- 4267 
st_crs(QUPR2.range)  <- 4267 
QUMI_clipped = st_intersection(QUMI.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
QUPR2_clipped = st_intersection(QUPR2.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = QUMI_clipped, col = "red", linewidth = 2)+
  geom_sf(data = QUPR2_clipped, col = "red", linewidth = 2)+
  geom_point(data = QUMI.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/QUMI.range.pdf", width = 12, height = 8)

Subset for Quercus minima. No range map available.

QUMI2.occ = gbif %>%
  filter(species == "Quercus minima") %>%
  select(species,decimalLatitude,decimalLongitude)
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_point(data = QUMI2.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/QUMI2.range.pdf", width = 12, height = 8)

Subset for Quercus muehlenbergii

QUMU.occ = gbif %>%
  filter(species == "Quercus muehlenbergii") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

QUMU.range = st_read("../../USTreeAtlas/shp/quermueh/")
## Reading layer `quermueh' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/quermueh' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 32 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -105.2417 ymin: 24.78449 xmax: -72.23889 ymax: 44.69553
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(QUMU.range)  <- 4267 
QUMU_clipped = st_intersection(QUMU.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = QUMU_clipped, col = "red", linewidth = 2)+
  geom_point(data = QUMU.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/QUMU.range.pdf", width = 12, height = 8)

Subset for Quercus nigra

QUNI.occ = gbif %>%
  filter(species == "Quercus nigra") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

QUNI.range = st_read("../../USTreeAtlas/shp/quernigr/")
## Reading layer `quernigr' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/quernigr' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 45 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -97.70858 ymin: 27.65694 xmax: -74.65406 ymax: 39.40351
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(QUNI.range)  <- 4267 
QUNI_clipped = st_intersection(QUNI.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = QUNI_clipped, col = "red", linewidth = 2)+
  geom_point(data = QUNI.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/QUNI.range.pdf", width = 12, height = 8)

Subset for Quercus pagoda. No range map available.

QUPA5.occ = gbif %>%
  filter(species == "Quercus pagoda") %>%
  select(species,decimalLatitude,decimalLongitude)
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_point(data = QUPA5.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/QUPA5.range.pdf", width = 12, height = 8)

Subset for Quercus palustris

QUPA2.occ = gbif %>%
  filter(species == "Quercus palustris") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

QUPA2.range = st_read("../../USTreeAtlas/shp/querpalu/")
## Reading layer `querpalu' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/querpalu' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 31 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -96.838 ymin: 34.436 xmax: -71.4752 ymax: 43.71538
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(QUPA2.range)  <- 4267 
QUPA2_clipped = st_intersection(QUPA2.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = QUPA2_clipped, col = "red", linewidth = 2)+
  geom_point(data = QUPA2.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/QUPA2.range.pdf", width = 12, height = 8)

Subset for Quercus phellos

QUPH.occ = gbif %>%
  filter(species == "Quercus phellos") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

QUPH.range = st_read("../../USTreeAtlas/shp/querphel/")
## Reading layer `querphel' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/querphel' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 42 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -96.26867 ymin: 28.91285 xmax: -73.14029 ymax: 40.74458
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(QUPH.range)  <- 4267 
QUPH_clipped = st_intersection(QUPH.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = QUPH_clipped, col = "red", linewidth = 2)+
  geom_point(data = QUPH.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/QUPH.range.pdf", width = 12, height = 8)

Subset for Quercus rubra

QURU.occ = gbif %>%
  filter(species == "Quercus rubra") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

QURU.range = st_read("../../USTreeAtlas/shp/querrubr/")
## Reading layer `querrubr' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/querrubr' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 155 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -97.78265 ymin: 30.7561 xmax: -59.68972 ymax: 49.08564
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(QURU.range)  <- 4267 
QURU_clipped = st_intersection(QURU.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = QURU_clipped, col = "red", linewidth = 2)+
  geom_point(data = QURU.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/QURU.range.pdf", width = 12, height = 8)

Subset for Quercus shumardii

QUSH.occ = gbif %>%
  filter(species == "Quercus shumardii") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

QUSH.range = st_read("../../USTreeAtlas/shp/quershum/")
## Reading layer `quershum' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/quershum' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 47 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -100.2109 ymin: 28.89789 xmax: -76.45997 ymax: 42.35025
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(QUSH.range)  <- 4267 
QUSH_clipped = st_intersection(QUSH.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = QUSH_clipped, col = "red", linewidth = 2)+
  geom_point(data = QUSH.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/QUSH.range.pdf", width = 12, height = 8)

Subset for Quercus sinuata. No range map available.

QUSIS.occ = gbif %>%
  filter(species == "Quercus sinuata") %>%
  select(species,decimalLatitude,decimalLongitude)
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_point(data = QUSIS.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/QUSIS.range.pdf", width = 12, height = 8)

Subset for Quercus stellata

QUST.occ = gbif %>%
  filter(species == "Quercus stellata") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

QUST.range = st_read("../../USTreeAtlas/shp/querstel/")
## Reading layer `querstel' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/querstel' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 74 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -101.3744 ymin: 27.61028 xmax: -69.95983 ymax: 41.76135
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(QUST.range)  <- 4267 
QUST_clipped = st_intersection(QUST.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = QUST_clipped, col = "red", linewidth = 2)+
  geom_point(data = QUST.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/QUST.range.pdf", width = 12, height = 8)

Subset for Quercus texana. No range map available.

QUTE.occ = gbif %>%
  filter(species == "Quercus texana") %>%
  select(species,decimalLatitude,decimalLongitude)
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_point(data = QUTE.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/QUTE.range.pdf", width = 12, height = 8)

Subset for Quercus velutina

QUVE.occ = gbif %>%
  filter(species == "Quercus velutina") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

QUVE.range = st_read("../../USTreeAtlas/shp/quervelu/")
## Reading layer `quervelu' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/quervelu' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 72 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -97.28247 ymin: 29.83254 xmax: -69.37499 ymax: 45.49608
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(QUVE.range)  <- 4267 
QUVE_clipped = st_intersection(QUVE.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = QUVE_clipped, col = "red", linewidth = 2)+
  geom_point(data = QUVE.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/QUVE.range.pdf", width = 12, height = 8)

Subset for Quercus virginiana

QUVI.occ = gbif %>%
  filter(species == "Quercus virginiana") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

QUVI.range = st_read("../../USTreeAtlas/shp/quervirg/")
## Reading layer `quervirg' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/quervirg' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 94 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -102.121 ymin: 22.41014 xmax: -75.71973 ymax: 37.46014
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(QUVI.range)  <- 4267 
QUVI_clipped = st_intersection(QUVI.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = QUVI_clipped, col = "red", linewidth = 2)+
  geom_point(data = QUVI.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/QUVI.range.pdf", width = 12, height = 8)

Subset for Robinia pseudoacacia

ROPS.occ = gbif %>%
  filter(species == "Robinia pseudoacacia") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

ROPS.range = st_read("../../USTreeAtlas/shp/robipseu/")
## Reading layer `robipseu' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/robipseu' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 12 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -95.69039 ymin: 32.40928 xmax: -76.87229 ymax: 41.12437
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(ROPS.range)  <- 4267 
ROPS_clipped = st_intersection(ROPS.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = ROPS_clipped, col = "red", linewidth = 2)+
  geom_point(data = ROPS.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/ROPS.range.pdf", width = 12, height = 8)

Subset for Salix caroliniana

CACA5.occ = gbif %>%
  filter(species == "Salix caroliniana") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

CACA5.range = st_read("../../USTreeAtlas/shp/salicaro/")
## Reading layer `salicaro' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/salicaro' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 100 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -100.292 ymin: 25.10583 xmax: -75.52333 ymax: 40.99319
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(CACA5.range)  <- 4267 
CACA5_clipped = st_intersection(CACA5.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = CACA5_clipped, col = "red", linewidth = 2)+
  geom_point(data = CACA5.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/CACA5.range.pdf", width = 12, height = 8)

Subset for Salix nigra

SANI.occ = gbif %>%
  filter(species == "Salix nigra") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

SANI.range = st_read("../../USTreeAtlas/shp/salinigr/")
## Reading layer `salinigr' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/salinigr' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 151 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -122.7956 ymin: 23.5848 xmax: -65.17285 ymax: 46.61042
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(SANI.range)  <- 4267 
SANI_clipped = st_intersection(SANI.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = SANI_clipped, col = "red", linewidth = 2)+
  geom_point(data = SANI.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/SANI.range.pdf", width = 12, height = 8)

Subset for Sassafras albidum

SAAL5.occ = gbif %>%
  filter(species == "Sassafras albidum") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

SAAL5.range = st_read("../../USTreeAtlas/shp/sassalbi/")
## Reading layer `sassalbi' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/sassalbi' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 74 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -96.55696 ymin: 28.49213 xmax: -69.92722 ymax: 44.81874
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(SAAL5.range)  <- 4267 
SAAL5_clipped = st_intersection(SAAL5.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = SAAL5_clipped, col = "red", linewidth = 2)+
  geom_point(data = SAAL5.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/SAAL5.range.pdf", width = 12, height = 8)

Subset for Sorbus americana

SOAM3.occ = gbif %>%
  filter(species == "Sorbus americana") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

SOAM3.range = st_read("../../USTreeAtlas/shp/sorbamer/")
## Reading layer `sorbamer' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/sorbamer' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 288 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -94.15333 ymin: 34.78817 xmax: -52.61445 ymax: 52.65512
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(SOAM3.range)  <- 4267 
SOAM3_clipped = st_intersection(SOAM3.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = SOAM3_clipped, col = "red", linewidth = 2)+
  geom_point(data = SOAM3.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/SOAM3.range.pdf", width = 12, height = 8)

Subset for Taxodium distichum (gbif) and Taxodium ascendens

TADI2.occ = gbif %>%
  filter(species == "Taxodium distichum") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

TADI2.range = st_read("../../USTreeAtlas/shp/taxodist/")
## Reading layer `taxodist' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/taxodist' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 85 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -99.82847 ymin: 24.54479 xmax: -75.04254 ymax: 38.92478
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(TADI2.range)  <- 4267 
TADI2_clipped = st_intersection(TADI2.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = TADI2_clipped, col = "red", linewidth = 2)+
  geom_point(data = TADI2.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/TADI2.range.pdf", width = 12, height = 8)

Subset for Thuja occidentalis

THOC2.occ = gbif %>%
  filter(species == "Thuja occidentalis") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

THOC2.range = st_read("../../USTreeAtlas/shp/thujocci/")
## Reading layer `thujocci' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/thujocci' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 134 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -102.1047 ymin: 35.59812 xmax: -63.26956 ymax: 55.0168
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(THOC2.range)  <- 4267 
THOC2_clipped = st_intersection(THOC2.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = THOC2_clipped, col = "red", linewidth = 2)+
  geom_point(data = THOC2.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/THOC2.range.pdf", width = 12, height = 8)

Subset for Tilia americana

TIAM.occ = gbif %>%
  filter(species == "Tilia americana") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

TIAM.range = st_read("../../USTreeAtlas/shp/tiliamer/")
## Reading layer `tiliamer' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/tiliamer' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 94 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -100.9682 ymin: 34.97041 xmax: -65.77465 ymax: 49.99133
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(TIAM.range)  <- 4267 
TIAM_clipped = st_intersection(TIAM.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = TIAM_clipped, col = "red", linewidth = 2)+
  geom_point(data = TIAM.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/TIAM.range.pdf", width = 12, height = 8)

Subset for Triadica sebifera. No range map available.

TRSE6.occ = gbif %>%
  filter(species == "Triadica sebifera") %>%
  select(species,decimalLatitude,decimalLongitude)
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_point(data = TRSE6.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/TRSE6.range.pdf", width = 12, height = 8)

Subset for Tsuga canadensis

TSCA.occ = gbif %>%
  filter(species == "Tsuga canadensis") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

TSCA.range = st_read("../../USTreeAtlas/shp/tsugcana/")
## Reading layer `tsugcana' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/tsugcana' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 180 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -93.74585 ymin: 33.5981 xmax: -59.68972 ymax: 48.33391
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(TSCA.range)  <- 4267 
TSCA_clipped = st_intersection(TSCA.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = TSCA_clipped, col = "red", linewidth = 2)+
  geom_point(data = TSCA.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/TSCA.range.pdf", width = 12, height = 8)

Subset for Tsuga caroliniana

TSCA2.occ = gbif %>%
  filter(species == "Tsuga caroliniana") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

TSCA2.range = st_read("../../USTreeAtlas/shp/tsugcaro/")
## Reading layer `tsugcaro' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/tsugcaro' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 9 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -83.89653 ymin: 34.66754 xmax: -79.44339 ymax: 37.66805
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(TSCA2.range)  <- 4267 
TSCA2_clipped = st_intersection(TSCA2.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = TSCA2_clipped, col = "red", linewidth = 2)+
  geom_point(data = TSCA2.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/TSCA2.range.pdf", width = 12, height = 8)

Subset for Ulmus alata

ULAL.occ = gbif %>%
  filter(species == "Ulmus alata") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

ULAL.range = st_read("../../USTreeAtlas/shp/ulmualat/")
## Reading layer `ulmualat' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/ulmualat' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 38 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -98.37525 ymin: 27.26967 xmax: -75.44749 ymax: 39.47826
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(ULAL.range)  <- 4267 
ULAL_clipped = st_intersection(ULAL.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = ULAL_clipped, col = "red", linewidth = 2)+
  geom_point(data = ULAL.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/ULAL.range.pdf", width = 12, height = 8)

Subset for Ulmus americana

ULAM.occ = gbif %>%
  filter(species == "Ulmus americana") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

ULAM.range = st_read("../../USTreeAtlas/shp/ulmuamer/")
## Reading layer `ulmuamer' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/ulmuamer' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 207 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -106.7759 ymin: 26.73745 xmax: -59.68972 ymax: 53.81448
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(ULAM.range)  <- 4267 
ULAM_clipped = st_intersection(ULAM.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = ULAM_clipped, col = "red", linewidth = 2)+
  geom_point(data = ULAM.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/ULAM.range.pdf", width = 12, height = 8)

Subset for Ulmus crassifolia

ULCR.occ = gbif %>%
  filter(species == "Ulmus crassifolia") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

ULCR.range = st_read("../../USTreeAtlas/shp/ulmucras/")
## Reading layer `ulmucras' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/ulmucras' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 15 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -100.5271 ymin: 25.47785 xmax: -89.93496 ymax: 36.42969
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(ULCR.range)  <- 4267 
ULCR_clipped = st_intersection(ULCR.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = ULCR_clipped, col = "red", linewidth = 2)+
  geom_point(data = ULCR.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/ULCR.range.pdf", width = 12, height = 8)

Subset for Ulmus rubra

ULRU.occ = gbif %>%
  filter(species == "Ulmus rubra") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

ULRU.range = st_read("../../USTreeAtlas/shp/ulmurubr/")
## Reading layer `ulmurubr' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/ulmurubr' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 50 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -101.2964 ymin: 29.64866 xmax: -68.63799 ymax: 48.01849
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(ULRU.range)  <- 4267 
ULRU_clipped = st_intersection(ULRU.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = ULRU_clipped, col = "red", linewidth = 2)+
  geom_point(data = ULRU.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/ULRU.range.pdf", width = 12, height = 8)

Subset for Ulmus thomasii

ULTH.occ = gbif %>%
  filter(species == "Ulmus thomasii") %>%
  select(species,decimalLatitude,decimalLongitude)

Read in range map

ULTH.range = st_read("../../USTreeAtlas/shp/ulmuthom/")
## Reading layer `ulmuthom' from data source 
##   `/Users/samanthaworthy/Documents/GitHub/USTreeAtlas/shp/ulmuthom' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 49 features and 5 fields
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: -99.79724 ymin: 35.06176 xmax: -72.06438 ymax: 47.69813
## CRS:           NA

The projection is geographic with latitude and longitude. The datum is NAD27 which is EPSG:4267. The ellipsoid is Clarke 1866, and the transformation parameters

st_crs(ULTH.range)  <- 4267 
ULTH_clipped = st_intersection(ULTH.range, states.map)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
ggplot()+
  geom_sf(data = states.map, fill = "white")+
  geom_sf(data = ULTH_clipped, col = "red", linewidth = 2)+
  geom_point(data = ULTH.occ, aes(x = decimalLongitude, y = decimalLatitude), color = "black", alpha = 0.7)+
  theme_classic()

ggsave("../Plots/ULTH.range.pdf", width = 12, height = 8)